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Abstract 

We consider the following problem: Given a rational matrix A E Q mxn and 
a rational polyhedron Q C M. m+P , decide if for all vectors b G W", for which 
there exists an integral z € V such that (b,z) £ Q, the system of linear in- 
equalities Ax ^ b has an integral solution. We show that there exists an 
algorithm that solves this problem in polynomial time if p and n are fixed. 
This extends a result of Kannan (1990) who established such an algorithm 
for the case when, in addition to p and n, the affine dimension of Q is fixed. 

As an application of this result, we describe an algorithm to find the maxi- 
mum difference between the optimum values of an integer program maxjcx : 
Ax ^ b, x £ Z"} and its linear programming relaxation over all right-hand 
sides b, for which the integer program is feasible. The algorithm is poly- 
nomial if n is fixed. This is an extension of a recent result of Ho§ten and 
Sturmfels (2003) who presented such an algorithm for integer programs in 
standard form. 

1 Introduction 

Central to this paper is the following parametric integer linear programming 
(PILP) problem: 

Given a rational matrix A £ Q mxn and a rational polyhedron Q C 
W n+P , decide if for all b £ W", for which there exists an integral 
z £ IP such that (b,z) £ Q, the system of linear inequalities Ax ^ b 
has an integral solution. 



1 



In other words, we need to check that for all vectors b in the set 



Q/Z p := {b G Q m : (b,z) G Q for some z G Z p } 



the corresponding integer linear programming problem Ax ^ b, x G Z n has a fea- 
sible solution. The set Q/Z p is called the integer projection of Q. Using this 
notation, we can reformulate PILP as the problem of testing the following V3- 
sentence: 

VbeQ/Z p 3xeZ n : Ax^b. (1) 

It is worth noticing that any polyhedron Q C M. m as well as the set of integral vec- 
tors in Q can be expressed by means of integer projections of polyhedra. Indeed, 

Q = Q/Z° and Qf]Z m = {(b,b) : b G Q}/Z m . 

In its general form, PILP be longs to the seco nd leve l of the polynomial hier- 
archy and is ITj -complete; see (|Stockmeyeiill976 ) and dWrathalj[l976l) . Kannan 



(|l990|) presented a polynomial algorithm to decide the sentence (OQ) in the case 
when n, p and the affine dimension of Q are fixed. This result was applied to de- 
duce a polynomial algorithm that solves the Frobenius problem when the number 
of input integers is fixed, see (|Kannanl . ll992|) . 

Kannan's algorithm proceeds in several steps. We informally describe it at this 
point as a way to decide V 3-statements (OQ) in the case p = 0. First Kannan pro- 
vides an algorithm which partitions the set of right-hand sides Q into polynomially 
many integer projections of partially open polyhedra Si,..., St, where each S ; is 
obtained from a higher-dimensional polyhedron by projecting out a fixed number 
of integer variables. Each 5/ is further equipped with a fixed number of mixed 
integer programs such that for each b G 5/ the system Ax ^ b is integer feasible, 
if and only if one of the fixed number of "candidate solutions" obtained from 
plugging b in these associated mixed integer programs, is a feasible integer point. 

To decide now whether (OQ) holds, one searches within the sets 5 ; individually 
for a vector b for which Ax ^ b has no integral solution. In other words, each of 
the candidate solutions associated to b must violate at least one of the inequalities 
in Ax ^ b. Since the number of candidate solutions is fixed, we can enumerate 
the choices to associate a violated inequality to each candidate solution. Each of 
these polynomially many choices yields now a mixed-integer program with a fixed 
number of integer variables. There exists ab G 5,- such that Ax ^ b has no integral 
solution if and only if one of these m ixed-integer pr ograms is feasible. The latter 
can be checked with the algorithm of lLenstral (|1983l) in polynomial time. 
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Contributions of this paper 



We modify the algorithm of Kannan to run in polynomial time under the assump- 
tion that only n and p are fixed. This is achieved via providing an algorithm 
that computes for a matrix A E Q mxn a set D C Z" of integral directions with 
the following property: for each b E W\ the lattice width (see Section [2]) of the 
polyhedron — {x : Ax ^ b} is equal to the width of this polyhedron along one 
of the directions in D. This algorithm is described in Section [3] and runs in poly- 
nomial time if n is fixed. The strengthening of Kannan's algorithm to decide 
V3-statements of the form © if n and p is fix ed follows then by using this result 
in the proof of Theorem 4.1 in jKannan , 1992h . 



We then apply this result to find the maximum integer programming gap for a 
family of integer programs. The integer programming gap of an integer program 

max{cx :Ax^b,xE Z"} (2) 

is the difference 

max{cx : Ax ^ b} — max{cx : Ax ^ b, x E Z"}. 

Given a rational matrix A E Q mxn and a rational objective vector c E Q'\ g(A,c) 
denotes the maximum integer programming gap of integer programs of the form 
©, where the maximum is taken over all vectors b, for which the integer pro- 
gram is feasible. Our algo rithm finds g(A,c) in polyno mial time if n is fixed. 



This extends a recent result of iHosten and Sturmfelsl (|2003|) . who proposed an al 



gorithm to find the maximum integer programming gap for a family of integer 
programs in standard form if n is fixed. 



Related work 



Kannan's algorithm is an extension of the polynomial al g orithm for integer lin 



Barvinok and Woods 



ear pr ogramming in fixed dimension by ILenstral (|1983|) 
d2003h presented an algorithm for counting integral points in the integer projec 
tion QjTLP of a poly tope Q C M. m+P . This algorithm runs in polynomial time if 
p and m are fixed, and uses Kannan's partitioning algorithm, which we extend 
in this paper. In particular, their algorithm can be applied to count the number 
of elements of the minimal Hilbert basis of a pointed cone in polynomial time 
if the dimension is fixed. We remark that a polyn omial test for the Hilbert basis 



property in fixed dimension was first presented by lCook et al.l (11984). Extensions 
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of Barvinok's algo rithm to compute counting function s for parametric polyhedra 
were pr esented in (IBarvinok and Pommersheiml 1 19991 : IVerdoolaege et all |2007|) 
and in (|Koppe and Verdoolaegel 2007 ). These counting functions are piecewise 
step-polynomials which involve roundup operations. With these functions at hand 
one can very efficiently compute the number of integer points in Pt, via evalua- 
tion at b. It is however not known how to use such piecewise step-polynomials to 
dec ide V3-statements effi c iently in fixed dimension. 



Hosten and Sturmfelsl (|2003l) proposed an algorithm to find the maximum in- 
teger programming gap for a family of integer programs in standard form, i.e., 
max {cx : Ax = b,x ^ 0,x £ Z"}. Their algorithm explo its short r a tional gen- 
erating functions for certain lattice point problems, cf. IBarvinok! dl994h and 
Barvinok and Woodsl ([2003), and runs in polynomial time if the number n of 



columns of A is fixed. However, the latter implies also a fixed number of rows 
in A, as we can always assume A to have full row rank. We would like to point out 
that our approach does not rely on rational generating functions at all. 



Basic definitions and notation 

For sets V and W in W and a number a we denote 

V + W := {v + w : v £ V, w £ W} and aW := {aw. w eW}. 

It is easy to see that if W is a convex set containing the origin and a ^ 1 , then 
aW CW. IfV consists of one vector v only, we write 

v + W := {v + w : w £ W} 

and say that v + W is the translate of W along the vector v. The symbol \a~\ de- 
notes the smallest integer greater than or equal to a, i.e., a rounded up. Similarly, 
[a\ stands for the largest integer not exceeding a, hence a rounded down. 

In this paper we establish a number of polynomial algorithms, i.e., algorithms 
whose running time is bounded by a polynomial in the input size. Following the 
standard agreements, we define the size of a rational number a = p/q, where 
p, q £ Z are relatively prime and q > 0, as the number of bits needed to write a in 
binary encoding: 

size(a) := 1 + \log(\p\ + 1)] + \log(q + 1)] . 

The size of a rational vector a = [a\, . . . ,a n ] is the sum of the sizes of its compo- 
nents: 

n 

size(a) := n + ^ size(a,-). 
(=l 
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At last, the size of a rational matrix A = [«,• 



G O rax " is 



size(A) :— mn + ^ ^2 

'=1.7=1 



size.(ajj) 



An open half-space in R" is the set of the form {x : ax < /3}, where a G R" is 
a row-vector and /3 is a number. Similarly, the set {x : ax ^ /3} is called a closed 
half-space. A partially open polyhedron P is the intersection of finitely many 
closed or open half-spaces. If P can be defined by means of closed half-spaces 
only, we say that it is a closed polyhedron, or simply a polyhedron. We need 
the notion of a partially open polyhedron to be able to partition the space (this is 
definitely impossible by means of closed polyhedra only). At last, we say that a 
partially open polyhedron is rational if it can be defined by the system of linear 
inequalities with rational coefficients and rational right-hand sides. 

Linear programming is about optimizing a linear function cx over a given 
polyhedron P in R w : 

max{cx : x G P} = — min{— cx : x G P}. 
If x is required to be integral, it is an integer linear programming problem 
max{cjc : x G PPlZ"} = -min{-cjc : x ePnZ n }. 



For details on linear and integer programming, we refer to ( Schrijver . 1986|) . Here 
we only mention that a linea r programming problem can be solved in polynomial 
time, cf. (|Khachiyanl . 1 1979Q . while integer linear programming is NP-complete. 
However, if the number of variables is fixe d, integer progr amming can also be 
solved in polynomial time, as was shown by iLenstral (|1983|) . Moreover, Lenstra 
presented an algorithm to solve mixed-integer programmi ng with a fixed num ber 
of integer variabl es. We remark that both algorithms — of lKhachiyan dl979h and 
of Lenstral ( 1983 ) — can be used to solve decision versions of integer and linear 
programming on partially open polyhedra. 

An integral square matrix U is called unimodular if |det(t/)| = 1. Clearly, if 
U is unimodular, then U~ l is also unimodular. A matrix of full row rank is said 
to be in Hermite normal form if it has the form [H 0] , where H — [/?,•_/] is a square 
non-singular non-negative upper-triangular matrix such that ha > hjj for all j > i. 
Given a matrix A of full row rank, we can find in polyno mial time a unimodu- 
lar m atrix U such that AU is in Hermite normal form; see (IKannan and Bacheml . 
1979Q . We remark that the Hermite normal form of an integral vector c is the 
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vector ae\, where a is the greatest common divisor of the components of c and 
e\ is the first unit vector. The unimodular matrix U such that cU = ae\ can be 
obtained directly while executing the Euclidean algorithm to compute the greatest 
common divisor. 



2 Flatness theorem 

We briefly review the algorithm to solve integer linear programming in fixed di- 
mension, as its basic ideas will be used in the following sections. Intuitively, if a 
polyhedron contains no integral point, then it must be "flat" along some integral 
direction. In order to make this precise, we introduce the notion of "lattice width." 
The width w c (K) of a closed convex set K along a direction c E R" is defined as 

w c {K) := max{oc : x E K} — min{oc : x E K}. (3) 

The lattice width w(K) of K (with respect to the standard lattice Z n ) is the mini- 
mum of its widths along all non-zero integral directions: 

w(K) := min{w c (^) :cEZ n \ {0}}. 

An integral row-vector c attaining the above minimum is called a width direction 
of the set K. Clearly, w(v + OtK) = ocw(K) for any rational vector v and any non- 
negative rational number a. Moreover, both sets K and v + CtK have the same 
width direction. 

Applications of the concept of lattice width in algorithmic number th eory and 



integ er programming rely upon the flatness theorem, which goes back to lKhinchin 



(1948) who first proved it for ellipsoids in W\ Here we state it for convex bodies, 



i.e., bounded closed convex sets of non-zero volume. 

Theorem 2.1 (Flatness theorem). There is a constant Co(n), depending only on n, 
such that any convex body K C W 2 with w(K) <X)(n) contains an integral point. 

The constant 6)(n) in Theorem l2.1l is referred to as the flatness constan t. The best 
known upper bound on C0(n) is 0(n 3 / 2 ), cf. ( Banaszczvk et al. L 1999b , altho ugh 



a linear dependence on n was conjectured, e.g., by Kannan and Lovasz jl988h . A 



linear lower bound on (d(n) was shown by lKantorl (LL999) and ISebd dl999) 



Throughout this paper we will mostly deal with rational polyhedra rather than 
general convex bodies. In this case, assumptions of non-zero volume and bound- 
edness can safely be removed from the theorem's statement. Indeed, if P C R" is a 
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rational polyhedron of zero volume, then it has width along an integral direction 
orthogonal to its (rational) affine hull. Further, let C be the characteristic cone of 
P: 

C = {y.x + y eP for allxEP}. 

If C = {0}, then P is already bounded. If C is full-dimensional, then the set x + C 
trivially contains an integral point, for any x G P (we can always allocate a unit 
box inside a full-dimensional cone). At last, if C is not full dimensional, then we 
can choose a sufficiently large box B CM" such that w(P) = w(P(~}B) and both 
P and P R B have the same width direction, which is orthogonal to the (rational) 
affine hull of C. If w(P) ^ (0(n), then PHB, and hence P, contains an integral 
point by Theorem 12. II 

How can we use this theorem to check whether a given rational polyhedron 
contains an integral point? The answer is in the following lemma, which is almost 
a direct consequence of the flatness theorem. 

Lemma 2.2. Let P C W 1 be a rational polyhedron of finite lattice width and let c 
be its width direction. Let 

j8 := min{cx : x e P} . (4) 
Then P contains an integral point if and only if the polyhedron 

PD{x:p^ cx^fi + (D(n)} 

contains an integral point. 

Proof. If w(P) < G)(n), then there is nothing to prove, since 

PCjjc:^ cx< /3 + «(/?)}. 

Suppose that w(P) ^ (o{n) and let P = y + Q, where y is an optimum solution of 
the linear program © and Q is the polyhedron containing the origin, 

Q:={x-y:xeP}. 

We denote 

Q':=$$Q and P>:=y + Q>. 

In other words, Q is P translated to contain the origin, Q' is obtained from Q by 
scaling it down, and P' is Q' translated back to the original position. It is easy to 
see that 

min{cx : x G P 1 } = cy = /3 . 
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cx = /3 cx = /3 + co{n) 



Figure 1: Illustration for the proof of Lemma [2721 

Since ^ 1 and Q is convex, we have Q' C Q. This implies P' C P. Yet, we 

have = w(Q), and therefore, w(P') = w(Q') = Co{n). 

By Theorem 12 .11 P' contains an integral point, say z. But then z also belongs 
to P and 

cz < max{cjc :iGP'} = /5 + 0)(n). 
This completes the proof. □ 

Suppose that we know a width direction c of a polyhedron 

P= {x:Ax^b} Ct". (5) 

Since c is integral, the scalar product cx must be an integer for any integral point 
x E P. Together with Lemma 12.21 it allows us to split the original problem into 
co{n) + 1 integer programming problems on lower-dimensional polyhedra 

Pn{x:cx= rjSl+7'}. 7 = 0,. 

where /3 is defined by ©. 

The components of c must be relatively prime, as otherwise we could scale 
c, obtaining a smaller width of P. Therefore its Hermite normal form is a unit 
row- vector e\. We can find a unimodular matrix U such that cU = e\, introduce 
new variables y := U~ l x and rewrite the original system of linear inequalities 
Ax ^ b in the form AUy ^ b. Since U is unimodular, the system Ax ^ b has an 
integral solution if and only if the system AUy ^ b has an integral solution. But 
the equation cx = |~j3] + j turns into eiy = |~j3] + j. Thus, the first component 
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of y can be eliminated. All together, we can proceed with a constant number of 
integer programming problems with a smaller number of variables. If n is fixed, 
this yields a polynomial algorithm. 

An attempt to generalize this approach for the case of varying b gives rise to 
the following problems. First, the width directions of the polyhedron © depend 
on b and therefore can also vary. Furthermore, even if a width direction c remains 
the same, it is not a trivial task to proceed recursively. The point is that /3, as it 
is defined in ©, also depends on b and the hyper-planes {x : cx = |~jB] + j} are 
not easy to construct with /3 being a function of b. In the following sections we 
basically resolve these two problems and adapt the above algorithm for the case 
of varying b. 

3 Lattice width of a parametric polyhedron 

A rational parametric polyhedron P defined by a matrix A E Q mXw is the family 
of polyhedra of the form 

P b := {x:Ax^b}, 

where the right-hand side b is allowed to vary over R' n . We restrict our attention 
only to those b, for which Pf, is non-empty. For each such b, there is a width 
direction c of the polyhedron Pj,. We aim to find a small set C of non-zero integral 
directions such that 

w{P b ) = min{w c (Pjb) : c E C} 

for all vectors b for which P b is non-empty. Further on, the elements of the set C 
are referred to as width directions of the parametric polyhedron P. It turns out that 
such a set can be computed in polynomial time when the number of columns in A 
is fixed. 

Let A E Q mxn be a matrix of full column rank. Given a subset of indices 

N = {h,...,i n } C {l,...,m}, 

we denote by An the matrix composed of the rows i\ , . . . , i n of A. We say that Af 
is a basis of A if An is non-singular. Clearly, any matrix of full column rank has 
at least one basis. Each basis Af defines a linear transformation 

F N :R m ^R'\ F N b=A N 1 b N , (6) 

which maps right-hand sides b to the corresponding basic solutions. We can view 
Fn as an n x m-matrix of rational numbers. If the point F^b satisfies the system 
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Ax ^ b, then it is a vertex of the polyhedron {x : Ax ^ b}. From linear program- 
ming duality we know that the optimum value of any feasible linear program 

max{cx : Ax ^ b} 

is finite if and only if there is a basis N such that c = yA^ for some row-vector 
y ^ 0. In other words, c must belong to the cone generated by the rows of matrix 
Apj. Moreover, if it is finite, there is a basis N such that the optimum value is 
attained at F^b. It gives us the following simple lemma. 

Lemma 3.1. Let P be a parametric polyhedron defined by a rational matrix A. If 
there exists a vector b' such that the polyhedron 

Py = {x : Ax ^ b'} 

has infinite lattice width, then w(Pt,) is infinite for all b. 

Proof. Suppose that the lattice width of P b is finite for some b and let c be a width 
direction. Then both linear programs 

max{cx : Ax ^ b} and min{oc : Ax ^ b} 

are bounded and therefore there are bases N\ and N2 of A such that c belongs to 
both cones 

Ci := {yA Nl : y > 0} and C 2 := {-yA Nl : y > 0} (7) 

generated by the rows of matrices A^ and — A^ 2 , respectively. But then the linear 
programs 

maxjoc : Ax ^ b'} and min{oc : Ax ^ b'} 
are also bounded, whence w c (Py) is finite. □ 

The above lemma shows that finite lattice width is a property of the matrix A. In 
particular Pq has finite lattice width if and only if P^ has finite lattice width for all 
b. Conversely, if Pq has infinite lattice width, then P^ also has infinite lattice width 
and therefore contains an integral point for all b. We can easily recognize whether 
Pq has infinite lattice width. For instance, we can enumerate all possible pairs of 
bases N\ and /V2 and check if the cones © have a common integral vector. Further 
we shall not deal with this "trivial" case and shall consider only those parametric 
polyhedra, for which w(Pq) is finite, and therefore w(Pb) is finite for any b. We 
say in this case that the parametric polyhedron P has finite lattice width. 
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Now, suppose that Pt, is non-empty and let c be its width direction. Then there 
are two bases N\ and N2 such that 

maxjcx: : Ax < b} = cF Nl b and min{cx : Ax ^ b} = cF Nl b (8) 

and c belongs to the cones C\ and C 2 defined by CO), see Figure[2] In fact, equations 
([8]) hold for any vector c in C\ fl C 2 . Thus, the lattice width of Pj, is equal to the 
optimum value of the following optimization problem: 

min{c(F M -F Nl )b : c e C\ nC 2 nZ n \{0}}. (9) 

The latter can be viewed as an integer programming problem. Indeed, the cones 
C\ and C2 can be represented by some systems of linear inequalities, say cD\ ^ 
and cD 2 ^ 0, respectively, where D\,D 2 E Z nxn . The minimum © is taken over 
all integral vectors c satisfying cD\ ^ and cD 2 ^ 0, except the origin. Since both 
cones C\ and C2 are simplicial, i.e., generated by n linearly independent vectors, 
the origin is a vertex of C\ fl C2 and therefore can be cut off by a single inequality, 
for example, cD\\ ^ —1, where 1 denotes the n-dimensional all-one vector. It 
is important that all other integral vectors c in C\ fl C2 satisfy this inequality and 
therefore remain feasible. Thus, the problem © can be rewritten as 

min{c(F Nl -F Nl )b : cD x sC 0, cD 2 ^ 0, cD { l < - lc G Z n }. 

For a given b, this is an integer programming problem. Therefore, the optimum 
value of © is attained at some vertex of the integer hull of the underlying poly- 
hedron 

{c:cDi ^0, cD 2 ^0, cD x \ ^ -1} (10) 
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Shevchenkol (|198l|) and lHayes and Larmanl (|1983|) proved that the number of ver- 
tices of the integer hull of a rational polyhed ron is polynomial in fix ed dimension . 
Tight bounds for number were presented in (ICook et al.l . ll992|) and (|Barany et al. . 
19921) . This gives rise to the next lemma. 



Lemma 3.2. There is an algorithm that takes as input a rational matrix A G Q mx " 
of full column rank, which defines a parametric polyhedron P of finite lattice 
width, and computes a set of triples (Fj,Gi,Ci) of rational linear transformations 
Fi, Gi : M, m — > W and a non-zero integral row-vector C{ G Z n (i = 1, . . . ,t) satisfying 
the following properties. For all b,for which Py is non-empty, 



(a) Fi and Gi provide, respectively, an upper and lower bound on the value of the 
linear function C[X in P^, i.e., for all i, 

CiGib ^ min{c,-* : x G Pb} < max{c,x : x G P^} < c,i^, 

(b) the lattice width ofPt, is attained along the direction cifor some i G { 1 , . . . , t} 
and can be expressed as 

w(P b ) = min Ci(Fj - Gi)b. 



(c) The number t of the triples satisfies the bound 

t^2m 2n (2n+l) n (24n 5 (t)) n -\ (11) 
where (j) is the maximum size of a column in A. 
The algorithm runs in polynomial time ifn is fixed. 

Proof. In the first step of the algorithm we enumerate all possible bases of A. 
Observe that there is at least one basis, since A is of full column rank. On the 
other hand, the total number of possible bases is at most m". Hence, the number 
of possible pairs of bases is bounded by m 2n . The algorithm iterates over all 
unordered pairs of bases and for each such pair {iVi,A^} does the following. 

Let C\ and C2 be the corresponding simplicial cones, defined by ©. These 
cones can be represented by systems of linear inequalities, cD\ ^ and c£>2 ^ 
respectively, where D\,D2 G Z nxn and the size of each inequality is bounded by 
4n 2 0, see (|Schrijveii 1 198a Theorem 10.2). As the origin is a vertex of the cone 
Ci flC2, it can be cut off by a single inequality; for example, cD{l ^ —1, where 
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1 stands for the n-dimensional all-one vector. The size of the latter inequality is 
bounded by 4n 3 0. 

Thus, there are exactly 2n+l inequalities in (flOl) and the size of each inequal- 
ity is bounded by 4n 3 0. This implies that the n umber of vertices of the integer 



hull of CUB is at most 2(2n + l)"(24n 5 0) n_1 , cf. dCook et all 1 19921). an d they all 
can be computed in polynomial time if n is fixed, cf. (|Hartmannl Il989|) . The al- 
gorithm then outputs the triple (F^ v F^ 2 ,c) for each vertex c of the integer hull of 
(fTOl) . where and F^ 2 are the linear transformations defined by ©. Since there 
are at most m 2n unordered pairs of bases and, for each pair, the algorithm returns 
at most 2(2« + \) n (24n 5 (j)) n ~ l triples, the total number of triples satisfies (TTTT) . 
as required. Parts (ja|) and © of the theorem follow directly from our previous 
explanation. □ 

The bound (PTTj) can be rewritten as 

for fixed n. Clearly, the greatest common divisor of the components of any direc- 
tion c, obtained by the algorithm must be equal to 1, as otherwise it would not be 
a vertex of the integer hull of (flOl) . This implies, in particular, that the Hermite 
normal form of any of these vectors is just the first unit vector e\ G W. 

It is also worth mentioning that if (Fj, G;,c;) is a triple attaining the minimum 
in Part © of Lemma [3T2l then we have 

w(Pb) ^ maxjcpc : x G Pb} — minjcpc : x G Pb} ^ ciFib — ciGib = w(Pb). 

Hence, Part fla]), when applied to this triple, turns into 

minjcpc : x G Pb} = C[G[b and maxjcpc : x G Pb} = c\F[b. 

For our further purposes, it is more suitable to have a unique width direction 
for all polyhedra Pb with varying b. In fact, using Lemma [3T2l we can partition the 
set of the right-hand sides b into a number of partially open polyhedra such that 
the width direction remains the same for all b belonging to the same region of the 
partition. 

Theorem 3.3. Let P be a parametric polyhedron of finite lattice width, defined 
by a matrix A G Q mxn of full column rank. Let Q C M. m be a rational partially 
open polyhedron such that Pb is non-empty for all b G Q. We can compute — in 
polynomial time, if n is fixed — a partition of Q into a number of partially open 
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polyhedra <2i , • -,Qt and, for each i, find a triple (F{, G,, c/) of linear transforma- 
tions Fi, Gt : W n — > M" and a non-zero vector C{ G II 1 , such that 

minjcpc : x G P b } = CjGjb, maxjcpc : x G P b } = CiFfi, 

and 

w{P b ) = w Ci (P h ) = Ci(F - Gi)b for all beQ,. 
If denotes the maximum size of a column in A, then t = 0(m 2n (j) n ~ l ). 



Remark. The statement of Theorem 13.31 is very analogous to (|KannanL 1 1992 . 
Lemma 3.1). However, there are several crucial differences. First, the number 
of regions in the partition obtained by Kannan's algorithm is exponential in n and 
the affine dimension jo of the polyhedron Q. Our algorithm yields a partition that 
is exponential in n only, hence polynomial if n is fixed. Also our algorithm runs 
in polynomial time if n is fixed but jo may vary. At last, the width directions Cj 
obtained by the Kannan's algorithm satisfy, for any b G Qu 

either w c .(P b ) < 1 or w Ci (P b ) < 2w(P b ) 

In contrast, we compute the exact width direction for each region in the partition. 
While the exact computation of width directions does not help much from the 
algorithmic point of view, removing the restriction on the dimension of Q turns 
out to be the main step towards the claimed generalization of Kannan's algorithm. 

Proof of Theorem [O] First, we exploit the algorithm of Lemma \32\ to obtain the 
triples (Fi,Gi,Ci), i = 1,...,?, with t = 0(m n (j) n ~ l ). These triples provide the 
possible width directions of the parametric polyhedron P. For each i = 1, . . . J, 
we define a partially open polyhedron Qi by the inequalities 

Ci(Fj-Gi)b <Cj(Fj-Gj)b, j= 1, 
diFi - Gi)b ^ Cj(Fj - Gj)b, ./ = /+l,...,f. 



Thus, 



min cj(Fj — Gj)b = c,(F ; - — G\)b 

i 



for all b G Qi. We claim that the intersections of the partially open polyhedra Qi 
with Q give the required partition. 

Indeed, let b G Q and let ji be the minimum value of Ci(Fi — Gj)b, i = l,...,t. 
Let / denote the set of indices i with Ci(Fi — Gi)b — fl. Then b G <2,- () , where j'q 
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is the smallest index in /. Yet, suppose that b G Q belongs to two partially open 
polyhedra, say Qj and Qj. Without loss of generality, we can assume i < j. But 
then we have 

Ci (Fi - Gt)b ^ Cj (Fj - Gj)b < a(Fi - Gi)b, 

where the first inequality is due to the fact b G Qj and the second inequality follows 
from b G Qf, both together are a contradiction. 

For the width directions, Lemma [X2l implies that 

w{P b ) = min Cj(Fj - Gj)b = c f (Fj - G { )b 
j 

for all b E Qifl Q. This completes the proof. □ 



Theorem 13.31 provides a unique width direction for each region Qj of the parti- 
tion. This resolves the first problem of adapting the algorithm for integer linear 
programming in fixed dimension to the case of varying b, which was addressed in 
the introduction. However, we still need to deal with the hyper-planes {x : cpc = 
|"j3] + j}, where /3 is the optimum value of the linear program minjcpc : x E Pb}. 
As mentioned above, /3 can be expressed as a linear transformation of b, namely 
/3 = ciGib. B ut |"j3] is no mo re a linear function of b, which makes the recursion 
complicated. Kannanl (1992) showed how to tackle this problem. We discuss it in 



the next section. 



4 Partitioning theorem and parametric integer pro- 
gramming 



The p roof of the following structural result follows from the proof of (IKannan , 



19921 Theorem 4.1) if it is combined with Theorem 13 .3[ 



Theorem 4.1. Let P be a parametric polyhedron of finite lattice width, defined by 
a rational matrix A G Q mxn of full column rank. Let Q C M m be a rational partially 
open polyhedron such that P b is non-empty for all b G Q. We can compute — 
in polynomial time, if n is fixed — a partition of Q into sets S\,...,St, and for 
each i, find a number of unimodular transformations Uij : W — > M" and affine 
transformations T\j : IR m — > IR'\ j = 1 , . . . , k t such that 

(a) each Si is the integer projection of a partially open polyhedron, Si = S'JlI'; 
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(b) for any b G 5,-, P^, fl Z n 7^ if and only ifP^, contains Uij \Tijb] for some index 

j; 

(c) the following bounds hold: 

t = 0((m 2n (j) n - 1 ) nls ^), h = 0(W(n)), fc ; = 0(2' l2 / 2 co(n)), i=l,...,f, 
where denotes the maximum size of a column in A and W(n) = YVi=\ G)(n). 

We do not repeat Kannan's proof here but give an intuition of why it is true in 
dimension 2. By Theorem [33] we can assume that the width-direction is invariant 
for all b E Q and by applying a unimodular transformation we can further assume 
that this width-direction is the first unit- vector e\, see Figure [3] 




x\ = \e\Gb\ x\ = \eiGb] +i 



Figure 3: Illustration of Theorem 14. II in dimension 2. 

Lemma [2T2l tells that Pf, contains an integral point if and only if there exists an 
integral point on the lines x\ = \eiGb~] + j for j = 0, . . . , a>(2). The intersections 
of these lines with the polyhedron P^ are 1 -dimensional polyhedra. Some of the 
constraints ax ^ /3 of Ax ^ b are "pointing upwards", i.e., aej < 0, where ej is 
the second unit-vector. Let a\x\ + a%xi ^ /3 be a constraint pointing upwards 
such that the intersection point (\e\Gb~\ + j,y) of a\xi -\-a2x2 = /3 with the line 
x\ = \eiGb~] + j has the largest second component. The line x\ = \eiGb~] + j 
contains an integral point in Pf, if and only if ( [eiG£>] + j, \y] ) is contained in P^. 
This point is illustrated in Figure [31 By choosing the highest constraint pointing 
upwards for each line x\ = \e\Gb \ + j, we partition the set of right-hand sides 
into polynomially many integer projections of partially open polyhedra. 
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In order to express the candidate solution (\e\Gb\ + j, \y]) in the form de- 
scribed in the theorem, observe that 



y = (j8 -a\x\)/a 2 . 

Since x\ is an integer and 

xi = \eiGb\ +j = e 1 Gb + j + y 

for some y G [0, 1), we can rewrite the equation \y] = |"(j8 — a\x\)/a 2 ~\ as 

[a\/a 2 \x\ + \y] = [j8/a 2 - {a\/a 2 }x i ] 

= \P/a 2 -{a\/a 2 }(e\Gb + j) - {a\/a 2 }f\, 

where {a\/a 2 } denotes the fractional part of a\/a 2 . Since {ai/a2}y lies between 
and 1, it suffices to check independently two different possibilities, namely, 

[ai/a 2 \x\ + \y] = \fi/a 2 - {a\/a 2 }{e\Gb + j)] 

and 

[ai/a 2 \xi + \y] = \P/a 2 - {a\/a 2 }(e\Gb + j) - 1] . 

Combined with x\ = \e\Gb~\ + j, each of the above equations yields a unimodular 
system with respect to the variables x\ and \y\ , with the right-han d side be i ng the 



round-up of an affine transformation of b. We refer the reader to (IKannanl 11992) 
to see the complete proof for arbitrary dimension. 

V 3 -statements 

Theorem 14 . 1 1 gives rise to a polynomial algorithm for testing sentences of the form 

VbeQ/Z p 3xeZ n : Ax^b, (12) 



when p and n are fixed. This algorithm was first described by IKannanl (| 1992) but 
he required, in addition, the affine dimension of Q to be fixed. Our improvement 
follows basically from the improvement in the partitioning theorem, while the 
algorithm itself remains exactly the same. We describe it here for the sake of 
completeness. First observe that we can assume that A has full column rank. 
Otherwise we can apply a unimodular transformation of A from the right to obtain 
a matrix [A' \ 0] , where A' has full column rank. 
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The idea is as follows: First we run the algorithm of Theorem 14.11 on input 
A and Q' C W n , where Q' is the set of vectors b, for which the system Ax ^ b 
has a solution. Then we consider each set Si returned by the algorithm of The- 
orem 14.11 independently. For each b £ 5 ; we have a fixed number of candidate 
solutions for the system Ax ^ b, defined via unimodular and affine transforma- 
tions as Uij\Tijb] . Each rounding operation can be expressed by introducing an 
integral vector: z = \Tjjb] is equivalent to Tjjb ^ z < Tip + 1. We need only a 
constant number of integer variables to express all candidate solutions plus a fixed 
number of integer variables to represent the integer projections 5, = S'J'L '. It re- 
mains to solve a number of mixed-integer programs, to which we also include the 
constraints (b,y) £ Q, y £ V . 

Theorem 4.2. There is an algorithm, that, given a rational matrix A £ Q mx " and 
a rational polyhedron Q C W m+P , decides the sentence (PT21) . The algorithm runs 
in polynomial time if p and n are fixed. 

Proof. Let P be a parametric polyhedron defined by the matrix A. First, we exploit 
the Fourier-Motzkin elimination procedure to construct the polyhedron Q' C M. m 
of the right-hand sides b, for which the system Ax ^ b has a (fractional) solu- 
tion. For each inequality ab ^ /3, defining the polyhedron Q', we can solve the 
following mixed-integer program 

ab>P, (b,y)eQ, yeZP, 

and if any of these problems has a feasible solution (y,b), then b is a vector in 
Q/Z p , for which the system Ax ^ b has no integral solution. Hence, we can 
terminate and output "no" (with b being a certificate). 

We can assume now that for all b £ QjTLP the system Ax ^ b has a fractional 
solution. By applying the algorithm of Theorem 14.11 we construct a partition of 
Q' into the sets Si,..., St, where each Sj is the integer projection of a partially 
open polyhedron, 5/ = S'JZ 1 ' . Since n is fixed, the /, are bounded by some con- 
stant, i = 1, . . . Furthermore, for each i, the algorithm constructs unimodular 
transformations t/,- ; - and affine transformations Ty, j = l,...,k(, such that Pj,, with 
b £ Su contains an integral point if and only if Utj \Tjjb] £ P\, for some j. Again, 
k\ is fixed for a fixed n, i= 1 , . . . , t. 

The algorithm will consider each index i independently. For a given i, 5,- can 
be described as the set of vectors b such that 
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has a solution for some integer z G Z'\ This can be expressed in terms of linear 
constraints, as S[ is a partially open polyhedron. Let xj = Uy \Tijb\ . The points xj 
can be described by linear inequalities as 

Tip C - / < Tj jh + 1 , xj = Ujjz j , 

where 1 is the all-one vector. Then Pi, does not contain an integral point if and only 
if Xj $l P\j for all j = 1 , . . . , . In this case, each Xj violates at least one constraint in 
the system Ax ^ b. We consider all possible tuples / of k[ constraints from Ax ^ b. 
Obviously, there are only m ' such tuples, that is, polynomially many in the input 
size. For each such tuple, we solve the mixed-integer program 

(b,y)eQ, (b,z)es' i: 

Tfjb < zj < Tijb + 1, j=l,...,ki, 

Xj = UijZ j, j = 1 , • • • , k{ , 

a ij x j > bij 7 j = 1 , . . . , kiy 

yeZP,zeZ l ', Zj eZ n , j = l,...,ki, 

where a{.x ^ b[. is the j-th constraint in the chosen tuple. Each such mixed-integer 
program can be solved in polynomial time since the number of integer variables 
is fixed (in fact, there are at most {ki +l)n + /,• integer variables). 

If there is a feasible solution b to one of these mixed-integer programs, then 
the answer to the original problem is "no" (with b being a certificate). If all these 
mixed-integer programs are infeasible, the answer is "yes". □ 

Rema rk. We wou ld li ke to point ou t that Theorem 14.21 can also be proved differ- 



ently. iBelll (|1977h andlScarfJ (| 19771) showed that if a system of linear inequalities 



Ax ^ b has no integral solution, then there is already a sub system of at m ost 2" 



inequalities that is infeasible in integer variables; see also (ISchrijveri 119861 The- 
orem 16.5). Applied to CQ) this means that there exists a b G QjV such that the 
system Ax ^ b has no integral solution, if and only there exist a b G Q/Z p and a 
subsystem A'x ^b'ofAx^b with at most 2" inequalities, which is infeasible in 
integer variables. Here b' is the projection of b into the according space. Since n 
is a constant, we can try out all Q") different subsystems of Ax ^ b and apply to 
each of these parametric polyhedra Kannan's algorithm to decide V 3- statements. 

However, a similar argument does not yield our extension (Theorem 14.11) of 
Kannan's partitioning theorem itself, which associates to each b a fixed set of 
candidate integer solutions, depending on the partially open polyhedron of the 
partitioning, in which b is contained. 
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5 Integer programming gaps 



Now, we describe how Theorem 14.21 can be applied to compute the maximum 
integer programming gap for a family of integer programs. Let A G Q rax " be a 
rational matrix and let c G Q n be a rational vector. Let us consider the integer 
programs of the form 

max{cx:Ax^b,xeZ n }, (13) 

where b is varying over R m . The corresponding linear programming relaxations 
are then 

maxjcjc : Ax ^ b}. (14) 
Consider the following system of inequalities: 

cx ^ /3, 
Ax^b. 

Given a vector b and a number /3, there exists a feasible fractional solution of 
the above system if and only if the linear program (fT4"l) is feasible and its value 
is at least /3. The set of pairs (/3,&) G M m+1 , for which the above system has 
a fractional solution, is a polyhedron in M. m and can be computed by means of 
Fourier-Motzkin elimination, in polynomial time if n is fixed. Let Q denote this 
polyhedron. 

Suppose that we suspect the maximum integer programming gap to be smaller 
than y. This means that, whenever /3 is an optimum value of (fT4l) . the integer 
program (fT3l must have a solution of value at least j8 — y. Equivalently, the system 

cO/3-y, (15) 
Ax ^ b, 

must have an integral solution. If there exists (£,/3) G Q such that (fT5l) has no 
integral solution, the integer programming gap is bigger than y. We also need to 
ensure that for a given b, the integer program is feasible, i.e., the system Ax ^ b 
has a solution in integer variables. 

Now, this is exactly the question for the algorithm of Theorem 14.21 Is there a 
(/3 , Z?) G Q' such that the system (fT5l) has no integral solution, but there exists y G 
Z" such that Ay ^ bl Here Q' = Q — y(l , 0) is the appropriate translate of the set <2. 
If the algorithm answers "no" with the certificate b, then the integer program (fT3~l) . 
with the right-hand side b, has no solution of value greater than /3 — y while being 
feasible. On the other hand, (/3 , Z?) G Q, thus the corresponding linear solution 
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has optimum value of at least /3. We can conclude that the maximum integer 
programming gap is greater than 7. This gives us the following theorem. 

Theorem 5.1. There is an algorithm that, given a rational matrix A G R mx ", a 
rational row-vector c G Q n and a number 7 checks whether the maximum integer 
programming gap for the integer programs (fT3l) defined by A and c is bigger than 
7. The algorithm runs in polynomial time if the rank of A is fixed. 

Using binary search, we can also find the minimum possible value of 7, hence the 
maximum integer programming gap. 
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